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Spacecraft reaction wheel maneuvers are limited by the maximum torque and/or 
angular momentum that the wheels can provide. For an w-wheel configuration, the torque or 
momentum envelope can be obtained by projecting the w-dimensional hypercube, 
representing the domain boundary of individual wheel torques or momenta, into three- 
dimensional space via the 3 xn matrix of wheel axes. In this paper, the properties of the 
projected hypercube are discussed, and algorithms are proposed for determining this 
maximal torque or momentum envelope for general wheel configurations. Practical 
strategies for distributing a prescribed torque or momentum among the n wheels are 
presented, with special emphasis on configurations of four, five, and six wheels. 


I. Introduction 

Many spacecraft employ more than three wheels for attitude control, both for redundancy and for the additional 
maneuvering or momentum storage capability of the extra wheels. The Swift Gamma Ray Burst Explorer, for 
example, employs six reaction wheels for rapid slewing to enable observation of the initial stages of ephemeral 
gamma ray burst events [1], and the James Webb Space Telescope (JWST) will employ a similar reaction wheel 
configuration [2]. In order to employ these wheels effectively, it is necessary to characterize the maximum torque 
and momentum available in any direction from the reaction wheels. This knowledge is also important in designing 
momentum unloading strategies. For brevity of exposition, we will emphasize angular momentum in the following. 
Fortunately, the mathematical analysis required to identify the maximum torque capability is identical to that 
required for angular momentum. This analysis is the subject of this paper. Previous work on the optimal use of more 
than three reaction wheels has largely been limited to the four-wheel case [3, 4]. This paper is an extension of an 
earlier conference presentation [5]. 


"Aerospace Engineer, Attitude Control Systems Engineering Branch, Code 591, AIAA Fellow 
f Senior GNC Engineer, Senior Member AIAA. 

* Principal Staff Engineer; GN&C. 
s Senior Scientist, Senior Member AIAA. 



In the most common case of nominally identical wheels with equal capabilities for both rotation directions, the 
available angular momenta fill the interior of an //-dimensional hypercube centered on the origin in reaction wheel 
angular momentum space. If the angular momentum capabilities vary from wheel to wheel, the angular momenta fill 
the interior of an //-dimensional rectangular parallelepiped in this space, where each side has length equal to twice 
the magnitude of the maximum angular momentum that a given wheel can provide. We will refer to this figure as a 
hypercube for brevity, at the expense of precision, since the case of equal wheel capabilities is the most common in 
practice. 

The available angular momenta in three-dimensional space fill a polyhedron given by the mapping of this 
hypercube by the 3 xn reaction wheel torque distribution matrix, whose columns are the unit vectors parallel to the 
reaction wheel axes in three-dimensional space. We will only consider “non-defective” reaction wheel 
configurations, by which we mean that no three of the reaction wheel spin axis directions are coplanar. As shown 
below, this polyhedron is the convex hull of the mapping of the vertices of the //-dimensional hypercube, which are 
the points at which ail n wheels are saturated in the positive or negative direction. The exterior vertices of the three- 
dimensional envelope are connected by edges that are the boundaries of the two-dimensional facets that form the 
envelope. 

We have developed algorithms for defining the envelope and for finding the maximum angular momentum in a 
given direction. These algorithms are also adapted to determine a distribution of an arbitrary angular momentum 
among the wheels that minimizes the maximum of the wheel momentum magnitudes. In this minimax distribution 
n - 2 wheels have the same maximum angular momentum magnitude, which means that the momenta lie on a scaled 
copy of the enveloping polyhedron whose size is determined by the required maximum wheel angular momentum 
magnitude rather than by the maximum wheel capability. The distribution algorithm incorporates a simple inner 
product test to determine which facet of the envelope gives the minimax magnitude, and hence which n ~ 2 wheels 
have the maximum magnitude. It is then straightforward to find the momenta of all the wheels necessary to produce 
the desired net angular momentum. We consider the specific examples of configurations of three, four, five, and six 
reaction wheels. 


II. Geometry of the Envelope 

The angular momenta of an awheel system can be represented by an /7-dimensional vector 

H (1) 

These angular momenta fill the interior of an //-dimensional hypercube, where the length of the side along the ith 
axis is two times the maximum angular momentum magnitude // . of the wheel along that axis. This can be the 

absolute maximum that the wheel can provide, or it can be some desired limiting momentum. The vertices of the 
hypercube are points where all wheels are saturated, i.e., they supply either their positive or negative maximum 
angular momenta. Each edge has all but one wheel saturated and is parallel to the axis of the unsaturated wheel. The 
facets of the hypercube have all but two wheels saturated. 

Denoting the wheel spin axis unit vectors by {w l w ? • • • w /7 } , the angular momentum vector in three- 
dimensional space is 

H = WH wheels (2) 

where 

W = [w, w, wj (3) 

is the angular momentum distribution matrix. The available angular momenta in three-dimensional space fill a 
polyhedron specified by the mapping of the /7-dimensional hypercube by Eq. (2). The vertices, edges, and facets of 
the polyhedron in three-dimensional space are mappings of the vertices, edges, and facets of the hypercube. 
However, this mapping is not 1:1, since some of the vertices, edges, and facets of the hypercube will map to the 
interior of the polyhedron in three-dimensional space. Thus the vertices of the available polyhedron in three- 
dimensional space are points where all wheels are saturated. The edges of the polyhedron have all but one wheel 



saturated, and an edge with wheel i unsaturated is parallel to that wheel’s spin axis w f - and has length 2 // . . The 

facets have all but two wheels saturated, and a facet with wheels / and j unsaturated is a parallelogram with sides 
parallel to w ; and and with normal vector in the direction of the cross product w / x w.. 


We can identify the bounding planes of the polyhedron, and thus determine the envelope, by considering all 
pairs i,j of wheels in turn. There are n(n-l)/2 such pairs for an //-wheel configuration. The /,y pair defines a set of 
2 th "' 2 facets of the hypercube in the /7-dimensional reaction wheel command space. Each facet corresponds to a 
saturated command for the n~~2 wheels other than wheels i and y, and with the angular momentum commands for 
wheels i and j varying over their full range over the facet. Under the transformation into three-dimensional space 
specified by Eq. (2), these facets map onto 2 n ~ 2 parallel planes with normal vector 

n = (w x w )/ w x w . (4) 

u v i J 7 1 ' J I v > 

On each of these planes, the wheel angular momentum in the direction of the normal to the plane is 

H “j, = i H k^k A y)= 2 CT *tfmax*(Vn^), ( 5 ) 

where o k = ±1 specifies the direction of saturation of the angular momentum on wheel k. Wheels i and j can clearly 
be omitted from the sum, since they don’t contribute to the angular momentum perpendicular to the plane. This is 
indicated algebraically by the vanishing of the triple product vv^ • (w ; . x w .) for k = iorj. T his triple product does 
not vanish for any other A, since we have assumed a “non-defective” reaction wheel configuration, for which no 
three of the reaction wheel spin axis directions are coplanar. 

The planes with the maximum and minimum angular momentum in the direction of the normal are the planes 
for which 

a k = sign [Wj, • (w x w y )] = sign (w* • n, ; ) (6a) 

or 

o k = -signfw* -(w. X w,)] = -sign^ • n tj ) (6b) 

for all k other than i and j. It is clear that these two planes and only these of the 2" 2 parallel planes for the i,j pair 
form part of the angular momentum envelope. On each of these bounding facets Hi and Hj vary over the range 
between ±H msxi and ±H m av , respectively. Saturation of H t and Hj give the four corners of a parallelogram, so we can 
identity eight of the vertices of the envelope for each i,j pair, corresponding to the two signs each for saturated H t 
and Hj , and the overall sign of the other saturated wheel angular momenta given by Eq. (6a) or (6b). 

Carrying out the above procedure over all pairs /, j will produce a total of An(n~\) vertices. These aren’t all 
distinct, however, for each vertex will be identified once for each bounding plane of the envelope that intersects 
there. This is useful information about the “structure” of the envelope, where “structure” denotes the number of 
vertices of the envelope, and the number of planes and edges that intersect at each vertex, as well as the specific 
information about which wheels are saturated in which direction at each vertex. It is interesting that this structure 
depends only on the unit vectors through their triple products w* -(w, xw y), and not on the maximum angular 
momentum capabilities of the wheels. If this angular momentum capability is changed, the shape of the envelope 
will change, but not the structure as defined here. It is obvious that the envelope has an inversion symmetry; for each 
point on the envelope with angular momentum H there is a corresponding point with angular momentum -H 
obtained by reversing all the wheel momenta. 

We will assume for the remainder of the paper that all the H maxi are equal, so we can omit the subscript / on 
H max . We also avoid the use of Eq. (6b) by distinguishing between the ij facet and the ji facet and always using Eq. 
(6a). Then on the ij facet of the polyhedron the wheel momentum is 


H: 


H l w { + ffw. 




max £ W^sign(W A 

k ij 




H,*, + HjWj 


H v . 

max ij ■ 


( 7 ) 



where 


v, y = E w i sign(w i .-n i;/ .). (8) 

k 5* ij 

The wheel momentum in the direction of the normal to the facet is 


where 


H ij = = A> = H ^ d i P 


^ = v -n = E V",, • 

fc^ij 


( 9 ) 

( 10 ) 


If the vector H f . n f . ends at a point on the hounding polyhedron, then ZT is the minimum angular momentum 

capability on that facet. This is not guaranteed to be the case, however, as we will show at the end of Section III.B of 
this paper. It is always true, though, that the minimum value of H } . over all the facets is the minimum angular 

momentum capability of the reaction wheel configuration. 


III. Optimizing the Angular Momentum Commands 


A. Optimality Criterion 


Optimal use of the angular momentum capability of a wheel configuration requires that equal demands be made 
on all wheels. In the case of wheels with identical capacities, this is accomplished by minimizing the maximum 
value, or L* norm, of the set of individual wheel angular momenta that will produce the desired net angular 
momentum. This minimax value H {) , which must be determined, determines the angular momentum envelope, as in 
the previous section, and the required angular momentum commands; it clearly cannot exceed the wheel capability 
of any wheel. The generalization to wheel assemblies with varying capabilities is straightforward, but we will not 
consider this uncommon case. 


B. Numerical Algorithm For Finding The Momentum Distribution 

Finding the L w momentum distribution requires identifying the facet of the wheel envelope that contains a point 
along the desired angular momentum axis. Assuming that we have found this to be facet ij , we write Eq. (7) as 

H = [w w ; \ ij ][H i Hj H 0 f, ( 11 ) 


where H , the maximum allowed momentum, has been replaced by H Q , the maximum momentum for this 
particular value of H. Equation (11) can be solved to give the wheel momenta as 


H, = [Vy ■ (w, x Wj )] ‘[(w J x Vy ) • H] 
^ — [Vy • ( W ( X Wy )] _l [( Vy x w ( ) • H] 
// 0 = WyH, 


where 


w = d n . 

V ij ij 


Note that the vectors v j; and are not unit vectors, and that v =— and w /( 
angular momenta for k ^ i, j are 


(12a) 

(12b) 

(12c) 

(13) 


Wy . The individual wheel 


H k = H 0 sign (w A ■ w (y ) for k * i,j . 


(14) 



The equations relating H wheeis to H are linear, so the mapping for facet ij can equivalently be written in matrix 
form as 


= H 0 5 ) 

Each matrix W* is a right inverse of W, meaning that their product is the 3x3 identity matrix; WW* = / 3x3 . These 

matrices can be pre-computed, and inspection of Eqs. (12)— (14) reveals that W*= W* , so it is only necessary to 
compute and store one matrix for each pair of parallel facets in this formulation. 


If spacecraft computer memory limitations preclude storing the W.. matrices, computations equivalent to Eqs. 
(12) — (14) must be performed onboard. We can avoid cross-product computations by using the identity 
I a 

|w, X W,| V.. = [w, ■ V.. -(w, • w y )(w. • v, ; )]w, + [w y • V.. -(w, ■ w y Xw, ■ V..)]w y + [(w ( x w y )-v.](w, X W j) (16) 
to rewrite Eqs. (12a) and (12b) as 


H, = \™i x j [w / ■ H - (w. * w. Xw y . • H)] - g y H 0 
i r 2 

Hj = |w ( . x w . | [ w j H - (w, • w j )(w, ■ H)] + g Jt H 0 

where 


(17a) 

(17b) 


I i-2 ^ a I a i~l a 

g = w x w [w. ■ v.. — (w. • w )(w * v. )] = w, x w n *(v x w .) . (18) 

\ i j\ 1 i tj v t j ij'* | i j j u y ij ,r v J 

r A i—2 

Storing the vectors w ;/ and the scalars wxwJ and for all the facets onboard along with the wheel axis unit 

vectors reduces the wheel angular momentum computations to only dot products, multiplications, and additions. 

The ij facet is a parallelogram with sides parallel to w ( - and w y . This is a subset of an infinite plane, defined by 

the conditions that |// ; .|< H 0 and |// ; | < H 0 . The point on this extended plane closest to the origin is the point 

where an angular momentum normal to the plane would lie. We want to find out whether this point lies on the facet, 
which we noted below Eq. (10) to not always be the case. Equations (17a) and (17b) show that H } = —g..H Q and 

Hj — gji^Q if H is normal to the plane. Thus the normal point lies on the facet if and only if |g^.| < 1 and |g /f | < 1 . 


C. Finding the Correct Facet 

l’f the computations of Eq. (12) yield a solution with either H f or H . having magnitude greater than H Q , it 

means that the ij facet is not the facet that gives the minimax momentum distribution. We now demonstrate a simple 
dot product test for finding the correct facet before performing these computations. Suppose the momentum vector 
H points to facet f p as shown in Figure 1. The normal vector nr. to the facet f }j and the momentum vector H, and 

thus the triangle OAG, lie in the plane of the figure. Consider any other facet fy whose normal is less than 90° from 
H. The normal vector n k} to this facet is not in the plane of the figure in general, and thus only the edge OF of the 

triangle OFE is guaranteed to be in the plane of the figure. The lengths \OA\ and |0i?| are given by 

\OA\ d (19a) 

|0£| = n..‘H, (19b) 

using Eq. (9) for H ip the distance of facet fj from the origin. The convexity of the momentum polyhedron means that 
| OF\ > \OG\ and therefore 

V H \OD\ \OC\JOC\ 1 OB | yH a yH 
^ax ff u \OE\ ! OF | | OG | \OA\ H tj H mm ' 



Equation (13) has been used for the outer terms in this inequality, which show that the correct facet is the one with 
the maximum value of w f . • H . The equality of W* and W* means that we only have to identify the correct pair 
of parallel facets, and the relation w . f =~~w // tells us that this is equivalent to finding the maximum value of 

|v H l- 


D. L t Alternative 


For comparative purposes, we consider the L 2 alternative to the L x analysis emphasized in this paper. The l 2 
norm minimizes the sum of the squares of the wheel momenta, rather than the maximum wheel momentum. It is 
well known that the wheel momentum distribution minimizing the L 2 norm is given by [6] 


H 


wheels 


W+ H, 


( 21 ) 


where W + is the Moore-Penrose pseudoinverse of W, which is given in this case by 

W + =W T (WW T y l . (22) 

Equation (21) will almost always result in the momentum of at of least one wheel having a greater magnitude than 
would result from Eq. (15). 


IV. Specific Cases 


We will consider four specific cases. The first three are symmetric configurations of three, four, and six wheels, 
respectively, with the wheel spin axes all at an angle t] from theyz plane, and uniformly distributed in azimuth about 
thex axis. The fourth case we will consider is the asymmetrical five-wheel case resulting from deletion of one wheel 
from the six-wheel configuration. The momentum envelopes for these four configurations are shown for comparison 
in Figures 2-5. The unit of length in these figures is the maximum angular momentum capability of a single wheel. 
The torque envelope, scaled to the maximum torque capability of a single wheel, is identical. Since the figures 
represent both envelopes, they are labeled Envelope, rather than Torque Envelope or Momentum Envelope. They all 
have the same elevation angle i] = 35.26° , which was chosen because it results in the largest minimum dimension of 
each of the envelopes, i.e. closest to spherical shape, as we will show algebraically below. For purpose of 
comparison, the four figures are all drawn to the same scale, viewed from the same vantage point, and have the same 
illumination sources. In a practical case for a specific mission, the choice of the angle tj is based on some 
compromise among convenience of fabrication, expected momentum accumulation, slewing requirements, and 
spacecraft inertia properties. 

A. 3-Wheel Configuration 


First consider the easy case of three wheels, with 


s 


W « 


0 


c 




(23) 


where s = si nrj and c ^ cos rj. In this case, there is 
momenta are given uniquely by the relation 

^whei 

which of necessity minimizes both norms. The structure of the envelope is that of a cube, with three planes meeting 
at each vertex. Thus each vertex will be identified three times by the procedure described above, and the number of 
distinct vertices is 4’3*(3-l)/3 - 8, the correct number of vertices of a cube. 


no distinction between the L and L, norms. The wheel 


i 





The maximum storage capacity on the 12 facet is found by computing 

Sc 


W l X W 2 


3s 

V3 s 


and 


d v =3cs/yj\ + 3s 2 


( 25 ) 


(26) 


The storage on ail the other facets is the same in this simple case. The storage as a function of r} is maximized for 

j] - sin” 1 (t/x/3 j - 35.26° , (27) 

which gives d n = 1 . For this value of ??, the axes are all perpendicular, so the momentum or torque envelope has the 

form of a tilted cube, as illustrated in Figure 2. It is not surprising that the minimum storage capacity for orthogonal 
wheel axes is the capacity of one wheel. It would be more common practice to orient three wheel axes along the 
spacecraft axes, in which case the wheel distribution matrix W would be the 3x3 identity matrix and the cube edges 
would be parallel to the jc, y, and z axes. The magnitude of the minimum momentum or torque capability does not 
depend on the orientation of the cube. 

B. 4-Wheel Configuration 

We next consider the often-employed case of four reaction wheels oriented in a tetrahedral configuration around 
the x axis. The angular momentum distribution matrix is 

s s s s 

W= 0 c 0 -c , (28) 

c 0 -c 0 

This four-wheel configuration has an envelope with 14 vertices. Six of the vertices are defined by the intersection of 
four facets each and 8 by three facets each. The torque or momentum envelope for this configuration, again 
assuming identical wheels, is shown in Figure 3. This envelope has a fourfold rotational symmetry about the x axis, 
which is the vertical axis in the figure. 

We now want to determine the storage capacity of the four-wheel configuration. In this case, there are only two 
different types of facet. Consider first the 12 facet. 


W 1 X W 2 = c 


so 


12 


= 4 cs/ 4 


1 + s 2 


(29) 


(30) 


The coefficients d i} for facets 23, 34, 41, 21, 32, 43, and 14 are the same. 


Now consider the 13 facet. 


so 


w ] x w 3 = 2 cs 


rfj 3 = 2c . 


(31) 


(32) 



The coefficients d { , for facets 24, 31, and 42 are the same. The storage as a function of is maximized at the angle 
that gives d ]2 = d u , which is 77 = 35.26°, the same value that maximizes storage in the three-wheel case. For this 

value of r/, the ratio of the minimum momentum storage to the capacity of a single wheel is 

d ]2 = d n = <Js/3 = 1.633 . 

The g . coefficients defined by Eq. (18) are 

§12 ~ § 23 = £34 ~ <§41 “ “§?] = ~§32 ” “§43 ~ — §14 “ (1 “ 3 )/(l +S ) (33a) 

<§14 “ §25 = §36 = §41 ~§52 = §63 = ^ ’ (33b) 

None of the g lJ coefficients in this four-wheel configuration has magnitude greater than unity for any wheel spin 
axis elevation angle, so every the facet contains its normal point. 


C. 6-Wheel Configuration 


The next case we consider has the wheels oriented with hexagonal symmetry around the x axis. The angular 
momentum distribution matrix is 




s s s 

o 4 c 4 c 

2 2 


c -c 
2 


-c -c 




s 


(34) 


This six-wheel configuration has an envelope with 32 vertices. Two of the vertices are defined by the intersection of 
six facets each, 18 by four facets each, and 12 by three facets each. Figure 4 illustrates the torque or momentum 
envelope for six wheels with identical capabilities in this configuration. This envelope has a sixfold rotational 
symmetry about the x axis. This is the reaction wheel configuration on the Swift or JWST spacecraft, except that the 
symmetry axis of the JWST wheel configuration is the z axis rather than thex axis. 

We now want to look at the special six-wheel instance of the general vector relations derived above. As we have 
seen, there are only three different types of facet. Consider first the 12 facet. 


Wj x w 2 


~y/3c 

s 

Ss 


so 


12 


: 6^3 cs/ >1: 


3 + r 


(35) 


(36) 


The coefficients for facets 23, 34, 45, 56, 61, 21, 32, 43, 54, 65, and 16 are the same. 


Now consider the 13 facet. 


so 


Wj x w 3 


c 

2 


Sc 

3 s 

4is 


13 



The coefficients d,, for facets 24, 35, 46,51, 62, 31 , 42, 53, 64, 15, and 26 are the same. 


(37) 


( 38 ) 



We finally consider the 14 facet. 


0 

1 

0 


(39) 


W, X w„ 


: 2 cs 


so 

d u ^l4ic . (40) 

The coefficients d {j for facets 25, 36, 41, 52, and 63 are the same. 

The minimum storage can be maximized by choosing 77 to maximize d 13 . This maximum is at 77 = 35.26° in 
this case, as in the three and four-wheel cases, which gives = b/ >/s — 2.683 , d ]3 = 8/3 = 2.667 , and 

d u = 2a/ 2 — 2.828 . The JWST angle 77 = 30° gives d n — 2.4962, d {3 — 2.6186 , and d u = 3, providing about 
7% less capacity in the worst direction than is provided by 77 = 35.26°. 


The g i} coefficients defined by Eq. (18) are 

£l2 =?23 =&34 =£45 =£56 = £61 = ~^21 =-£32 =-£43 =-£54 =-£65 “ "*16 = ( 3 “ 1 1 ^ )/( 3 + ^) (41a) 

£ 13 = £ 2 4 = £ 3 5 - £46 = £51 - £62 " -£i3 = “£42 = -£53 “ -£ 6 4 = -£15 “ -£26 =( 3 - 7 r)/(l + 3 s 2 ) (41 b) 

£14 “ £25 = £36 = £41 =£52 = £63 = 0 • ( 4 1 c ) 


Because all the g tj coefficients of Eq. (41c) have magnitude less than unity, the six facets described by that 
equation all contain their normal points for any wheel elevation angle. The twelve facets described by Eq. (41a) 
contain their normal points only if s 2 <,3/5 . The twelve facets described by Eq. (41b) contain their normal points 

only if s 2 s> l/ 5 . In the range l/ 5 £ s 2 £ 3/5 , or 26.565° ^ ?/ £ 50.768° , all the facets contain their normal points and 
all the vectors H n fJ end on the bounding polyhedron. 

D. 5 -Wheel Configuration 

We will not consider a symmetric configuration of five wheels, but rather the asymmetric configuration 
resulting from the failure of one wheel from the symmetric six-wheel configuration. Figure 5 shows the torque or 
momentum envelope for the case where the angular momentum distribution matrix is obtained by deleting one 
column from the matrix of Eq. (34). 

There are six different types of facets in the failed wheel cases. We will only consider the case that wheel 6 has 
failed; the other cases can be found easily by permutation of indices. In this case 


d n = d 7] = d i5 = d 5i = 5-y/J cs/^3 + s 2 

(42a) 

d ]3 ~ d 3X = d 3 5 * d~ 3 = 6 cs / Vl + 3s 2 

(42b) 

d \4 “ ^41 = d 25 « d 5 2 = 3^3 cj 2 

(42c) 

d 23 = d n = d 34 = d 43 = 4 a/5" cs/ / l + s 2 

(42 d) 

d l5 = d 5] - lcs/ 4 1 + 3s 2 

(42e) 

d 2A “ d^ 2 ** Scs/4l + 3 ^ 2 

(42i) 


For the wheel inclination angle that gives the maximum momentum storage in the worst direction, 77 = 35.26°, the 
smallest of these is d 24 = d A2 , and it has just 5/8 of the smallest value for the full six-wheel configuration. 



The g.j coefficients defined by Eq. (18) satisfy the equations 


g 2 \ =-£45 “2-&I2 = 2 + Ssa =105 2 /(3+^ 2 ) 

(43a) 

£32 = ~&34 =1_ ^23 =1+ &I3 = %S 1 /0 + S 1 ) 

(43b) 

§\1 ~ ~§5i ~ 1 — £31 “ 1 + £35 — 3c / (1 + 35") 

(43c) 

g 24 =-g 42 =2(l-2s 2 )/(l + 3s 2 ) 

(43d) 

g 5] =-g :5 =2(\-4s 2 )/(l + 3s 2 ) 

(43 e) 

225 =-£41 = 1+ g 52 =1 ~Sl4 - ] /4. 

(43 f) 


There is much less symmetry among these coefficients in the failed wheel case. Examination of Eq. (43a) shows that 
in this configuration, all of the vectors H fJ n tj can end on the bounding polyhedron only if the spin axis elevation 

angle has the value rj « 35.26° , and the other equations show that this condition is sufficient. 

V, Avoiding Equal Wheel Speeds 

It is often undesirable to have multiple wheels running at the same speed, so in the case of n > 3 wheels we can 
use vectors in the («-3)-dimensional null space of W to separate the wheels speeds. Angular momentum commands 
in the null space can also be used to drive a wheel’s speed rapidly through zero or past the frequency of a flexible 
mode of the spacecraft structure without perturbing the pointing of the spacecraft. 

The null space is one-dimensional in four-wheel case, so all null vectors are multiples of 

p = in - 1 i-if. (44) 

Even though there is only one null vector, the Z, x algorithm assigns momenta on every facet with signs that enable 
this null vector to separate any equal wheel speeds. 

Several convenient null vectors in the six-wheel case are: 


u 0 =-£[l-l 1-1 1-1]' 

(45a) 

u 14 =i[0 1-1 0 I -If 

(45 b) 

c 

N> 

U1 

il 

wh 

T 

0 

0 

(45c) 

u 36 -i[l-l 0 1-1 Of 

(45d) 

u.=I[0 1-3 4-3 if 

(45e) 

U =I[1 0 1-3 4-yf 

1 6 

(451) 

u 3 = I [—3 10 1-3 4 f 

(45g) 

u 4 =i[4-3 1 0 l-3f 

(45h) 

u 5 =|[-3 4-3 1 0 if 

(451) 

u 6 ={[1-3 4-3 1 Of 

(45j) 


The indices indicate the locations of zeros in the vectors. The normalization factors are chosen to give unit 
vectors, but they’re really irrelevant because there’s no requirement for the null space vectors to be normalized. Note 
that u 14 4- u 25 -f u 36 = 0 and that u Q is orthogonal to u ]4 , u ?5 , and u 36 . 



In the six-wheel case, u Q and any two of u 14 , u 25 , or u 36 constitute a basis for the three-dimensional null 
space. The failed wheel cases are equally simple. In the case of a failed wheel 6, for example, u 6 and u 36 constitute 
an orthogonal basis for the two-dimensional null space. The other failed wheel cases are similar. 

An alternative method of avoiding wheels running at the same speed is to distribute the momenta according to a 
blend of the L m and L 2 distribution laws, 

= + ( 46 ) 

where A; is a scalar constant between 0 and 1. Choosing the smallest value of k that gives an acceptable separation of 
wheel speeds maximizes the momentum storage capacity. This alternative is really a special case of using vectors in 

the null space to separate the wheel speeds, because all the columns of W* -W + lie in the null space. 


VI. Slews 


A slew can be executed under reaction wheel control by varying the total reaction wheel momentum vector in 
three-dimensional space according to some slew law specifying a torque H . The L distribution of this torque 
among the reaction wheels is given in analogy with Eq. (15) by 


H 


wheels 


W/H, 


(47) 


where the ij facet is the one with the maximum value of w y . • H . Note that this torque will not necessarily keep the 

reaction wheels in a state that minimizes the L m norm of their momenta, because the facet on which * H is 

maximized is not the same as the facet on which w ;/ • H is maximized, in general. This can be dealt with by setting 

up reaction wheel control loops to feedforward the L x torque and then to use tachometer feedback to asymptotically 
drive the wheels to the L m momentum distribution. These loops can be incorporated into a wheel drag compensation 
loop as illustrated schematically in Figure 6. 


Another option is to maintain the Loo momentum distribution by commanding the wheel torques to be the 
derivatives of the momenta computed by Eq. (15), This protects against momentum saturation of the wheels and is 
easier to implement, but it gives discontinuous wheel torque commands when the momentum moves from one facet 
of the polyhedron to another, and it will generally result in larger wheel torque commands than those computed by 
Eq. (47). Spacecraft requirements differ, and the desire to minimize wheel momentum must be balanced against the 
desire to minimize wheel torques for any given application. Generally speaking, the time required to execute small 
slews tends to be dominated by the torque capability of the reaction wheels, and the time for large slews by the 
momentum capacity of the wheels. 


The L 2 distribution law does not have this complication; taking the derivative Eq. (21) gives 


H 


wheels 


W + H 


(48) 


because the pseudoinverse is a constant matrix. Thus applying torques distributed according to an L 2 law will 
preserve an L 2 momentum distribution. 


VII. Conclusions 

We have presented algorithms for determining the envelope of available torques and angular momenta in three- 
dimensional space for general configurations of n reaction wheels. These fill a polyhedron given by projection into 
three-dimensional space of an ^-dimensional hypercube in reaction wheel space. This polyhedron is the convex hull 
of the projection of the vertices of the hypercube, the points at which all the wheels are saturated in the positive or 
negative direction. We have also presented an algorithm for finding the optimal angular momentum commands to 
produce any desired angular momentum. The optimality criterion is that these commands minimize the maximum 
value, or L* norm, of the vector of individual wheel angular momenta. The algorithm also determines the maximum 



angular momentum available in any given direction. The exact solution algorithm is computationally feasible for 
onboard processing, especially if the necessary quantities for a specific wheel configuration are precomputed. The 
onboard computations could be further reduced by storing the maps of vector angular momentum or torque to 
individual wheels as a set of nx 3 matrices, one for each pair of opposite facets, at the cost of greater onboard 
computer storage. 
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